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The signatures from Sunyaev-Zeldovich effects and free-free emission in the intergalactic and intracluster 
medium and at galactic scales probe the structure evolution at various cosmic times. The detection of these 
sources and, possibly, the precise imaging of their spatial structure at the high resolution and sensitivity achiev- 
able with SKA will greatly contribute to the comprehension of crucial cosmological and astrophysical aspects, as 
the physical conditions of early ionized halos, quasars and proto-galactic gas and the mapping of the thermal and 
density structure of clusters of galaxies, with an extremely accurate control of extragalactic radio sources, which 
represents also a useful by-product for future CMB spectrum experiments devoted to probe the thermal plasma 
history at early times. 



1. INTRODUCTION 

The thermal plasma in the intergalactic and in- 
tracluster medium and at galactic scales leaves 
imprints on the cosmic microwave background 
(CMB) through the Thomson scattering of CMB 
photons on hot electrons and the free-free emis- 
sion. Although not specifically devoted to CMB 
studies, because of its high resolution and the lim- 
ited high frequency coverage, the extreme sensi- 
tivity and resolution of SKA may be fruitfully 
used for a detailed mapping of the above effects 
on dedicated sky areas. A particularly accurate 
control of the extragalactic radio sources repre- 
sents a useful by-product for future CMB spec- 
trum experiments devoted to probe the thermal 
plasma history at high redshifts. The SKA high- 
est frequencies (~ 20 GHz) are, of course, the 
most advantageous to detect Sunyaev-Zeldovich 
effects on the CMB because of the steeper de- 
crease of the synchrotron radio emission with the 
frequency. 



2. THERMAL AND KINETIC 

SUNYAEV-ZELDOVICH EFFECT 
TOWARDS GALAXY CLUSTERS 

The scattering of CMB photons from hot elec- 
trons in galaxies and clusters of galaxies produces 
a frequency dependent change in the CMB bright- 
ness. If the hot electron gas is globally at rest 
with respect to the observer only the thermal 
Sunyaev-Zeldovich (SZ) effect [ED (see also [12]) 
will be present; differently, a bulk peculiar mo- 
tion, V r , of the hot electron gas produces a kinetic 
SZ effect. In the Rayleigh- Jeans (RJ) region the 
first effect produces a decrement of the surface 
brightness, AI t h, towards the cluster. The sec- 
ond effect produces either a decrement or an in- 
crement, Alk, depending on the direction of the 
cluster velocity with respect to the observer. Ne- 
glecting relativistic corrections: 

AI th = I yg(x) (1) 

and 

AI k = -I (V r /c)T e h(x) , (2) 
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where J = {2h P /c 2 ){k B T CMB /h) 3 . Here 

T e = y n e a T dl (3) 

and 

y = J (k B T e /m e c 2 )n e a T dl (4) 

are respectively the Thomson optical depth and 
the Comptonization parameter [^Ej integrated 
over the cluster along the line of sight, n e being 
the electron density, and 

h{x) = x 4 e x /(e x - If (5) 

g(x) = h(x)[x(e* + l)/(e x -l)-4}, (6) 

where x — hpv jk B TcMB is a dimcnsionless pho- 
ton frequency. The two effects have a differ- 
ent frequency dependence that in principle allows 
their separation through multi-frequency observa- 
tions. SKA observations in the RJ regime (where 
h(x) r~j g(x) — > x 2 ) can be combined with milli- 
metric observations (in particular g(x) ~ and 
h(x) is maximum at ~ 217 GHz). 

Many thousands of galaxy clusters will be iden- 
tified by XMM (~ 10 3 ), Planck (~ 10 4 ), and 
SDSS (~ 510 5 ). The typical angular sizes of 
galaxy clusters range from ~ arcmin to few tens 
of arcmin. 

The SKA sensitivity and resolution mainly de- 
pends on the used array collecting area. By con- 
sidering a frequency band of ~ 4 GHz at 20 GHz, 
the whole instrument collecting area will allow 
to reach a (rms) sensitivity of ~ 40 nJy in one 
hour of integration with an angular resolution of 
~ 1 mas (considering ~ 3000 km maximum base- 
line). By using only about 50% of the collecting 
area within ~ 5 km, the (rms) sensitivity in one 
hour of integration is ~ 80 nJy with a resolution 
of ~ 0.6". 

Although the major role on the study of the SZ 
effects towards galaxy clusters will be played by 
dedicated telescopes operating at ~ arcmin res- 
olutions with frequency coverages up to ~ milli- 
metric wavelengths [23J > with the 50% of the SKA 
collecting area it will be possible to accurately 
map the SZ effect of each considered cluster, 
particularly at moderately high redshifts, with 



an extremely precise subtraction of discrete radio 
sources (see also Sect.0J. The combination with 
X-ray images, in particular with those expected 
by the wide field imager (WFI) on board the 
X-ray Evolving Universe Spectroscopy (XEUS) 
satellite 1 by ESA designed to reach a resolution 
of 0.25" on a FOV of 5' - 10', will allow to accu- 
rately map the thermal and density structure of 
the gas in galaxy clusters. 

Finally, we note that it is also possible to study 
the SZ effect (both thermal and kinetic) from 
clusters in a statistical sense, namely through 
its contribution to the angular power spectrum, 
Ce, (see Appendix A) of the CMB secondary 
anisotropies. This topic has been investigated in 
several papers (see, e.g., |35l54ll 514511 ll2l^ . At 
sub-arcmin scales (i.e. at multipoles t £ 10 4 ) sec- 
ondary anisotropies from thermal (more impor- 
tant at £ ,$ few x 10 4 ) and kinetic (more impor- 
tant at £ £ few x 10 4 ) SZ effect dominate over 
CMB primary anisotropy whose power signifi- 
cantly decreases at multipoles £ £ 10 3 because 
of photon diffusion (Silk damping effect |47p. 
Their angular power spectrum at I ~ f0 4 — 
10 5 (w f0~ 12 — f0~ 13 in terms of dimensionless 
C(.l{2£ + l)/4-7r) could be in principle investigated 
with the sensitivity achievable with SKA (see Fig. 
6 and Appendix A). On the other hand, at the 
SKA resolution and sensitivity the contribution 
to fluctuations from foreground sources (both dif- 
fuse radio emission, SZ effects, and free-free emit- 
ters) at galaxy scales probably dominates over the 
SZ effect from clusters. 

3. THERMAL SUNYAEV-ZELDOVICH 
EFFECT AT GALAXY SCALE 

The proto-galactic gas is expected to have a 
large thermal energy content, leading to a de- 
tectable SZ signal, both when the protogalaxy 
collapses with the gas shock-heated to the virial 
temperature |4fl55| . and in a later phase as the 
result of strong feedback from a flaring active nu- 
cleus (see, e.g., |1<)I30I31I1I38T2T| . The astrophys- 
ical implications of these scenarios have been re- 
cently investigated by 

1 http: / / sci.esa.int / science-e / www/area/index. cfm?fareaid=25 
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A fully ionized gas with a thermal energy den- 
sity e gas within the virial radius 
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transfers to the CMB an amount Ae ~ 
(e g asAc)2(i?vir/c) of thermal energy through 
Thomson scattering producing a Comptonization 
parameter j^E] y — (l/4)Ae/ecMB- Here h is the 
Hubble constant in units of 100 kms _1 Mpc _1 
and p vir ~ 200p„, p u = 1.88 x 10" 29 /i 2 (l + 
z v - lr ) 3 g cm -3 is the mean density of the uni- 
verse at the virialization redshift z v j r ; ccmb = 
obbTcmb - 4.2xl0 _13 (l+z) 4 ergcm" 3 , a B B be- 
ing the black-body constant and Tcmb = Tq(1 + 
z) the temperature of the CMB; t c is the gas cool- 
ing time by Thomson scattering. 

Assuming the binding energy (-Eb.gas = 

M gas W 2 ir , V vil = 162^ 1 /3( 1 + z )l/2( Mv . r/10 'l2 MQ) l/3 

km s _1 being the circular velocity of the galaxy at 
its virial radius |32l4j ) to characterize the thermal 
energy content of the gas, E gas , the amplitude of 
the SZ dip in the RJ region can be written as: 
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This SZ effect shows up on small (typically sub- 
arcmin) angular scales. 

Quasar-driven blast-waves could inject into the 
ISM an amount of energy several times higher 
than the gas binding energy, thus producing 
larger, if much rarer, SZ signals. A black-hole 
(BH) accreting a mass -Mbh with a mass to radi- 
ation conversion efficiency €bh releases an energy 
i?BH = cbh^bhc 2 . We adopt the standard value 
for the efficiency £bh = 0.1 and assume that a 
fraction fh = 0.1 of the energy is fed in kinetic 
form and generates strong shocks turning it into 
heat. 



Using the recent re-assessment by [21] of the 
well known correlation between the BH mass and 
the stellar velocity dispersion, Mbh = 1-4 x 
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The amplitude of the SZ dip in the RJ region due 
to quasar heating of the gas is then estimated as: 
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Following |38|. we adopt an isothermal density 
profile of the galaxy. The virial radius, encom- 
passing a mean density of 200p u , is then: 
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corresponding to an angular radius: 
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where (1a(z) is the angular diameter distance. 

The angular scales of these SZ signals from 
galaxies are of the order of « 10", then of particu- 
lar interest for a detailed mapping with the SKA 
and XEUS in the radio and X-ray, respectively. 
The probability of observing these SZ sources on 
a given sky field at a certain flux detection level 
and the corresponding fluctuations are mainly de- 
termined by the redshift dependent source num- 
ber density 4>sz(Ssz, z) per unit interval of the 
SZ (decrement) flux 5s z- The lifetime of the con- 
sidered SZ sources is crucial to determine their 
number density. 

For quasar-driven blast-waves the lifetime of 
the active phase, tgz> is approximately equal 
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to the time for the shock to reach the outer 
boundary of the host galaxy. Assuming a self- 
similar blast-wave expanding in a medium with 
an isothermal density profile, p oc r~ 2 , we have: 

*» - ^^rx^r (i3) 

The evolving B-band luminosity function of 
quasars, 4>(Lb, z), can be then adopted to derive 
the source number density 4>sz(Ssz, z) according 
to 

0sz(5 , sz,^) = 4>{L B ,z)-^- ^-5- , (14) 
tq.opt u&sz 

where Lb{Ssz,z) is the blue luminosity of a 
quasar at redshift z causing a (negative) SZ flux 
Sszj and £ q , pt is the duration of the optically 
bright phase of the quasar evolution. 

For the proto-galactic gas isz should be re- 
placed by the gas cooling time, t CO oi- Assum- 
ing that quasars can be used as effective sign- 
posts for massive spheroidal galaxies in their early 
evolutionary phases ^5] an d that they emit at 
the Eddington limit and using the relation by 
|12| between the mass of the dark-matter halo, 
M v i r , and the mass of the central black- hole, 
Mbh/10 8 M q ~ O.l(M vir /lO 12 M ) 1 - 65 , the num- 
ber density of sources with gas at virial tem- 
perature can be straightforwardly related to the 
quasar luminosity function <\>{L*b , z) . 

In spite of the many uncertainties of these 
models, it is remarkable that the CMB fluctua- 
tions (dominated at small scales by the Poisson 
contribution) induced by the SZ effect of these 
source populations could contribute to the CBI 
anisotropy measure and, in particular, could ex- 
plain the angular power spectrum (see also Ap- 
pendix A) found by BIMA at multipoles I w 
(4-10) x 10 3 (see Fig.HJ. 

A direct probe of these models and, possibly, 
their accurate knowledge through a precise high 
resolution imaging is then of particular interest. 
Fig. [21 shows the number counts at 20 GHz pre- 
dicted by these models: in a single SKA FOV 
about few x 10 2 — 10 3 SZ sources with fluxes above 
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Figure 1. Angular power spectrum of SZ effects 
at 30 GHz compared to CMB primary fluctuation 
power spectrum and CBI (2Z| (box) and BIMA 
[H] (data points) measures. Solid lines represent 
clustering (bottom line), Poisson (middle line) 
and global (upper line) contributions from quasar 
driven blast- waves. Dashed lines represent clus- 
tering (bottom line at high £), Poisson (middle 
line at high £) and global (upper line) contribu- 
tions from proto-galactic gas. The latter are actu- 
ally upper limits since, because of the uncertainty 
in the cooling time, the extreme assumption that 
icooi = texp has been adopted in the computation. 
Dots refer to the overall contribution. 



^100 nJy could be then observed in few hours of 
integration. Given the typical source sizes, we ex- 
pect a blend of sources in the SKA FOV at these 
sensitivity levels, while much shorter integration 
times, ~ sec, on many FOV would allow to obtain 
much larger maps with a significant smaller num- 
ber of resolved SZ sources per FOV. Both surveys 
on relatively wide sky areas and deep exposures 
on limited numbers of FOV are interesting and 
easily obtainable with SKA. 

A different scenario to jointly explain the power 
excess found by BIMA and the high redshift 
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Figure 2. Number count predictions at 20 GHz 
for SZ effects as function of the absolute 
value of the flux from proto-galactic gas heated 
at the virial temperature (dashes) assuming 
Mgas/Mvir = 0.1 and from quasar driven blast- 
waves (solid line). The exponential model for the 
evolving luminosity function of quasars is derived 
by |37| for an optical spectral index of quasars 
a = 0.5 (S v oc v a ). The parameters have been set 
at e B H = 0.1, f h = 0.1, fcboi = 10, £ q ,opt = 10 7 yr. 



reionization detected by WMAP (z reion ~ 
15 — 20) has been proposed by It involves 

hot gas winds powered by pair-instability super- 
novae (SN) explosions from the first generation of 
very massive stars at very low metallicity able to 
photoevaporate the gas in the halo potential. The 
SN remnants should then dissipate their energy in 
the intergalactic medium (IGM) and about 30- 
100% of their energy would be transferred to the 
CMB via Compton cooling. The resulting SZ 
effect from this sources is relevant in statistical 
sense. It is claimed to explain the high I BIMA 
excess of the CMB angular power spectrum and 
to be able to generate a global Comptonization 
distortion parameter y ~ few x 10~ 6 . However, 
it is estimated to be too faint (~ few x 10~ 2 nJy) 



to be observable even by SKA. 

4. SKA CONTRIBUTION TO FUTURE 
CMB SPECTRUM EXPERIMENTS 

The current limits on CMB spectral distor- 
tions and the constraints on energy dissipation 
processes |Ae/e,| 10 -4 in the plasma @3] are 

mainly set by COBE/FIRAS |2BIT3| . CMB spec- 
trum experiments from space, DIMES J23] (see 
also U) at A > 1 cm and FIRAS II [H] at 

A ;$ 1 cm, have been proposed with an accuracy 

potentially able to constrain (or probably detect) 
energy exchanges 10-100 times smaller than the 
FIRAS upper limits. In particular, experiments 
like DIMES may probe dissipation processes at 
early times (z }t 10 5 ) resulting in Bose-Einstein 

like distortions |50I8I5| and free-free distortions 
|2] possibly generated by heating (but, although 
disfavoured by WMAP, in principle also by cool- 
ing [471] ) mechanisms at late epochs (z ;$ 10 4 ), be- 
fore or after the recombination era The high 
redshift reionization detected by WMAP J2j sup- 
ports the existence of late coupled Comptoniza- 
tion (with Ae/ei ~ Ay few x 10~ 6 , e$ being 
the CMB energy density before the beginning of 
the dissipation process) and free-free distortions 
(with a highly model dependent amplitude) [5]. 
Typical distorted spectra potentially detectable 
by DIMES are shown in Fig. |3] To firmly ob- 
serve such small distortions the Galactic and ex- 
tragalactic foreground contribution should be ac- 
curately modelled and subtracted. 

Recent progress on radio source counts at 
1.4 GHz have been presented in PH]. On the 
other hand, the very faint tail of radio source 
counts is essentially unexplored and their con- 
tribution to the radio background at very low 
brightness temperature is not accurately known. 
For illustration, by assuming differential source 
number counts, N(S), given by logN(S)/AN ~ 
alogS + b, with AiV ~ 150S- 2 5 sr" 1 Jy" 1 (5 in 

Jy) [12, for a - 0.4 - 0.6 and b (0.5 - 1), 

we find a contribution to the radio background 
at 5 GHz from sources between ~ 1 nJy and 
~ 1/iJy between few tens of /iK and few mK. 
These signals are clearly negligible compared to 
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Figure 3. CMB distorted spectra as functions 
of the wavelength A (in cm) in the presence of 
a late energy injection with Ae/e^ = 5 x 10 -6 
plus an early/intermediate energy injection with 
Ae/ei = 5 x 10~ 6 occurring at the "time" Comp- 
tonization parameter yh = 5,1,0.01 (from the 
bottom to the top; in the figure the cases at 
Uh = 5 and 1 are indistiguishable at short wave- 
lengths; solid lines) and plus a free-free distortion 
with Ub = 10~ 6 (dashes), yh is defined by Eq. (4) 
with dl — cdt and T e = Tqmb when the integral 
is computed from the time of the energy injection 
to the current time. 



the accuracy of current CMB spectrum experi- 
ments, in particular a,t X > 1 cm, but are signifi- 
cant at the accuracy level on CMB distortion pa- 
rameters potentially achievable with experiments 
like DIMES. This effect is small compared to the 
Galactic radio emission, whose accurate knowl- 
edge currently represents the major astrophysical 
problem in CMB spectrum experiments, but, dif- 
ferently from Galactic emission, it is isotropic at 
the angular scales of few degrees and can not be 
then subtracted from the CMB monopole tem- 
perature on the basis of its angular correlation 
properties. With accurate absolute measures on a 



wide frequency coverage a fit including both CMB 
distorted spectra and astrophysical contributions 
can be searched (see 03] for an application to 
FIRAS data) but a direct radio background es- 
timate from precise number counts will certainly 
improve the robustness of this kind of analyses. 

The SKA sensitivity at 20 GHz will allow the 
detection (to 5c) of sources down to a flux level 
of ~ 200 nJy (~ 60, 20, 6 nJy) in 1 (10, 10 2 , 10 3 ) 
hour(s) of integration over the ~ 1 mas (FWHM) 
resolution element; similar numbers (from ~ 250 
to 8 nJy in an integration time from 1 to 10 3 
hours, respectively) but on a resolution element 
about 10 times larger will be reached at « GHz 
frequencies by using a frequency bandwidth of 
about 25%. 

Therefore, the SKA accurate determination of 
source number counts down to very faint fluxes 
can directly help the solution of one fundamental 
problem of the future generation of CMB spec- 
trum space experiments at A ^ 1 cm. 

5. FREE-FREE EMITTERS 

The epoch of reionization has been extensively 
studied in the recent years and in particular since 
the WMAP 1-yr data release. Reionization af- 
fects the CMB both in anisotropics at large and 
small scales and in the spectrum. The under- 
standing of the ionizing emissivity of collapsed 
objects and the degree of gas clumping is cru- 
cial for reionization models. The observation of 
diffuse gas and Population III objects in thermal 
bremstrahhlung as a direct probe of these quanti- 
ties has been investigated by j3J3 ■ Free- free emis- 
sion produces both global and localized spectral 
distortion of the CMB. A natural way to distin- 
guish between free-free distortion by ionized halos 
rather than by diffuse ionized IGM is represented 
by observations at high resolution of dedicated 
sky areas and by the fluctuations in the free-free 
background. In the model by [33] halos collapse 
and form a starburst lasting t = 10 7 yr, then 
recombine and no longer contribute to the free- 
free background. By adopting a Press-Schechter 
model |40I3| for the number density of collapsed 
halos per mass interval, dnps/dM, |33| exploited 
the expression by 0S] f° r the collapse rate of halos 
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per mass interval per unit comoving volume: 

dJ yform 1 dDdn PS/nr . SI 



dM 



a 2 (M)D- 



5(15) 



here D(z) is the growth factor and S c = 1.7 is 
the threshold above which mass fluctuations col- 
lapse. The expected comoving number density of 
ionized halos in a given flux interval as a function 
of redshift 



dN] 



halo 



dSdV 



(S,z) 



*( z )-*° dN? ™ 

dt T7T- 

dM 



dM 



t(z) 



dS 



(16) 



can be then computed given the expected flux 
from a halo of mass M at redshift z, S — S(M, z) 1 
and the starburst duration, t Q . Adopting a cut- 
off mass for a halo to be ionized of M* = 10 8 (1 + 
z/10)- 3 / 2 M Q (the critical mass needed to attain 
a virial temperature of 10 K to excite atomic 
hydrogen cooling), |33| computed the number 
counts of sources above the flux limit S c from 
the zeroth moment of the intensity distribution 
moments due to sources above a redshift z m j n , 
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,„ dN halo dV 
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by setting S max -> oo and S mm (z) = 
max.(S c , S*(z)), where S*(z) denotes the flux 
from a halo of minimum mass at redshift z. 

The relation iV rccomb = a B (n 2 )V w (1 - 
fesc)Ni on , between the production rate of re- 
combination line photons, iV rec0 mb, and the pro- 
duction rate of ionizing photons, N- lon , (here 
as is the recombination coefficient and f esc 
(w some%) is the escape fraction for ionizing 
photons) implies that the source luminosities in 
Ha and free- free emission (cx n 2 e V) are propor- 
tional to the production rate of ionizing pho- 
tons. Over a wide range of nebulosity conditions 
[T5| found that ~ 0.45 Ha photons are emitted 
per Lyman continuum photon; thus L(Ha) = 
1.4 x 10 41 7V ion /(10 53 phs" 1 )ergs" 1 . Given the 
free-free volume emissivity 03] in the case of 
an approximate mild temperature dependence 
with a power law (a velocity averaged Gaunt 
factor §ff — 4.7 is assumed), e v = 3.2 x 
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Assuming the starburst model of J7j normalised 
to the observed metallicity 10~ 3 Z© < Z < 
10 -2 i?Q of the IGM at z ~ 3 (resulting into a con- 
stant fraction of the gas mass turning into stars, 
1.7% < f s tar < 17%, in a starburst which fades 
after ~ 10 7 yr), [^j derived a production rate of 
ionizing photons as a function of halo mass given 
by: 



N ioD (M) = 2 x 10 
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which specifies the above free-free ionized halo 
luminosity. The corresponding flux is then: 
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Clearly, SKA will allow to detect only bright 
sources with deep exposures. The ionized halo 
number counts can be calculated from Eq. (17). 
The result by [S| is reported in Fig. H SKA 
should be able to detect ~ 10 4 individual free- free 
emission sources with z > 5 in 1 square degree 
above a source detection threshold of 70 nJy. The 
redshift information from the Balmer line emis- 
sion detectable by the Next Generation Space 
Telescope (NGST) can be used to discriminate 
ionized halos from other classes of radio sources. 

Ionized halos may contribute to the tempera- 
ture fluctuations. In particular, the Poisson con- 
tribution is predicted to be larger (smaller) than 
the clustering one at scales smaller (larger) than 
~ 30" [HSj- On the other hand, both are likely 
dominated by the radio source contribution. 

Finally, the integrated emission from ionized 
halos produces a global CMB spectral distortion, 
ATff — c 2 {S) /2kBV 2 , that can be computed 
from the mean sky averaged signal (S) . By using 
Eq. (17) (since no point source removal is feasi- 
ble at degree scales) with z min and S c = 0, [3*3*| 
found a free-free distortion AT// = 3.4 x 10~ 3 K 
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Figure 4. Number of sources which may be de- 
tected in the 1° by SKA, as a function of the 
threshold flux S c . Realistic limiting fluxes for 
point source detection are shown. The extrap- 
olated source counts from are also shown. 
From ES]. 



at 2 GHz, corresponding to a free-free distortion 
parameter ys — 1.5 x 10~ 6 , well within the ob- 
servational capability of DIMES [2SCj- 

6. CONCLUSION 

Although not specifically devoted to CMB 
studies, because of its high resolution and the 
limited high frequency coverage, the extreme sen- 
sitivity and resolution of SKA may be fruitfully 
used for a detailed mapping of the thermal plasma 
properties in the intergalactic and intracluster 
medium and at galaxy scale. 

On the basis of the existing literature, we dis- 
cussed the SKA contribution to the understand- 
ing of some topics of relevant interest in cosmol- 
ogy, from the detailed mapping of the SZ effect 
towards cluster of galaxies, to the number counts 
and imaging possibilities of the SZ effect from 
early quasars and proto-galactic gas and of the 



free-free emission from early ionized halos. We 
also discussed the contribution of SKA to the fu- 
ture CMB spectrum space experiments, devoted 
to the comprehension of the plasma thermal his- 
tory at early epochs, through the extremely accu- 
rate control of extragalactic radio source counts 
at very faint fluxes. 

A. SKA SENSITIVITY IN TERMS OF 
ANGULAR POWER SPECTRUM 

The properties of a large population of astro- 
physical objects can be studied by exploiting the 
details of the sources in dedicated fields and by 
analyzing the field statistical properties. 

The statistics of temperature anisotropy is typ- 
ically analyzed in spherical harmonics m : 

jr(l)=Y, E "imY tm W, (21) 

1=1 m=-£ 

where 7 is the observation unit direction vector. 
By assuming isotropy around the observer, ai m 
should have zero mean, (a^ m ) = 0, and variance 
Ce 2 given by: 

Oe = (\a 2 em \) = ^ n: J2 a im- ( 22 ) 

m 

The correlation function of the temperature 
anisotropy, C{9) = (^-(71) ^-(72)), is related to 
the angular power spectrum by the equation: 

C(9) = ^-y(2l + l)C e P e (cos9)- (23) 

47T — ' 

/ 

here cos 9 = 71 • 72 and is the Legendre poly- 
nomial. 

At SKA resolution and sensitivity the standard 
approaches to compute the angular power spec- 
trum of anisotropics from classes of astrophysical 

2 Roughly speaking, a multipolc ~ I corresponds to an 
angular scale i9/deg = 180/^. In this convention Ci is 
dimensionless. It is also useful to bear in mind that the 
same expansion in spherical harmonics can be considered 
for the (physical) temperature fluctuation ST, instead of 
the dimensionless temperature fluctuation ST/T; in this 
case the physical dimension of Ci will be the square of a 
temperature. Given the CMB monopole temperature of 
2.725 K the dimensionless C e will be ~ 7.4 X 10 12 

smaller than the C< expressed in terms of fiK 2 (in ther- 
modynamic temperature). 
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\10 2 FOV-, 



Figure 5. Relative uncertainty on CMB angular 
power spectrum recovery from combined cosmic 
and sampling variance for some representative sky 
coverages. 



ag m coefficients are random variables (possibly 
following a Gaussian distribution), at a given t, 
and therefore their variance, Cg, is \ 2 distributed 
with 2£+l degrees of freedom. The relative vari- 
ance SCi on Ci is equal to a/2/ {21 + 1) which 
is quite relevant at low I because of the rela- 
tively small number of available modes. This is 
the so-called "cosmic variance" that defines the 
ultimate limit on the accuracy at which a given 
model defined by an appropriate set of param- 
eters can be constrained by the angular power 
spectrum. Another similar variance in anisotropy 
experiments is related to the sky coverage since 
the detailed anisotropy statistical properties may 
depend on the considered sky patch. This vari- 
ance depends on the observed sky fraction, / s ky 
At the largest multipoles achievable with a given 
experiment the most relevant uncertainties are re- 
lated to the experiment resolution and sensitivity. 
All these terms contribute to the final uncertainty 
on the angular power spectrum according to |22| : 



/sky(2£+f) 



Aa 2 



NC e Wi 



(24) 



sources, typically based on the assumption that 
sources above the detection level do not signif- 
icantly merge in the same image pixel, are no 
longer valid since the expected number of (no 
longer point-like) sources per resolution element 
may be large and the source diffuse emission ex- 
tends over several (or many) image pixels (see, 
e.g., Fig. Further studies are then necessary to 
characterize the correlation properties of source 
populations as observed at SKA resolution and 
sensitivity. 

In spite of this uncertainty, it is important to 
characterize the nominal SKA capability to re- 
cover these kinds of statistical information. 

We then "rewrite" here the typical SKA sensi- 
tivity in terms of accuracy in the recovery of the 
angular power spectrum. 

Since each given anisotropy field is a single real- 
ization of a stochastic process, it may be different 
from the average over the ensemble of all possible 
realizations of the given (true) model with given 
parameters. This translates into the fact that the 



where A is the size of the surveyed area, a is 
the rms noise per pixel, N is the total number 
of observed pixel, and Wg is the beam window 
function. For a symmetric Gaussian beam Wi = 
exp(-^(£+ l)cr|) where a B = FWHM/V8hr2 de- 
fines the beam resolution. In the limit of an ex- 
periment with infinite sensitivity (er = 0) 5Ci/Ci 
is determined only by cosmic and sampling vari- 
ance; in this limit, by assuming also full sky cov- 
erage (/sky = 1), SCe/Cg is determined only by 
the cosmic variance. 

At ~ 20 GHz a typical SKA field of view (FOV) 
is of ~ 4' x 4' and a sky coverage of ~ I0 2 FOV 
corresponds to about ~ 40' x 40', possibly achiev- 
able with mosaicing techniques. 

The relative error SCgjCi due to cosmic and 
sampling variance is shown in Fig. [S] by consid- 
ering a binning of ~ 10% in I which allows to 
reduce SCg /Q by a factor y/0A£ and is not crit- 
ical, being quite smooth the Ce dependence on 
£ for the kinds of anisotropics considered in this 
work. As evident the sky coverage correspond- 
ing to - I0 2 SKA FOV (or to - I SKA FOV) 
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Figure 6. Absolute uncertainty on dimensionless 
CMB angular power spectrum recovery due to the 
instrumental noise for some reference cases (see 
also the text). 



allows to reach a relative "fundamental" uncer- 
tainty less than ~ 10% at £ ~ 10 4 (or at ~ 10 5 ) 
and significantly smaller at larger multipolcs. 

The smallest angular scales (or the highest mul- 
tipolcs) achievable with a good accuracy are set 
by the instrument resolution and sensitivity. 

By using about 50% of the collecting area 3 
SKA would allow the mapping of diffuse 
anisotropies up to £ ~ 10 6 . We will use this case 
as a reference. 

The absolute error bCg due to the instrumental 
noise is shown in Figs. El and \7\ again by consid- 
ering a binning of ~ 10% in £, respectively in 
terms of Cg£(2£ + 1)/4tt with dimensionless Cg 

3 In this context, the whole SKA collecting area could 
represent a real improvement only provided that surface 
brightness (or temperature) anisotropy maps at higher res- 
olution could be degraded at lower resolution to improve 
the rms (white noise) sensitivity per resolution element 
(according to the ratio of the corresponding FWHM) by 
averaging the brightness temperatures of the small pix- 
els which constitute a larger pixel, a point that should be 
accurately verified through detailed simulations. 



Figure 7. The same as in Fig. but in terms of 
the temperature fluctuation y / Cg£(2£ + 



and of y/Ci£(2£+l)/Air with C t in ^K 2 (note 
that we do not refer here to the noise angular 
power spectrum but to the uncertainty on the sig- 
nal angular power spectrum recovery, i.e. to the 
rms of "residual" noise angular power spectrum 
after the subtraction of its expectation value ob- 
tained on the basis of the accurate knowledge of 
the instrumental noise properties, as necessary in 
CMB anisotropy experiments). 



ACKNOWLEDGMENTS It is a pleasure 
to thank A. Cavaliere, L. Danese, C. Fanti, 
F. Finelli, G.L. Granato, A. Lapi, N. Man- 
dolesi, P. Platania, L.A. Popa, R. Salvaterra, and 
L. Silva for collaboration and discussions. The 
use of the CMBFAST code is acknowledged. 

REFERENCES 

1. Aghanim N., Balland C, Silk, J., 2000, A&A, 
357, 1 

2. Barlett J.G., Stebbins A., 1991, ApJ, 371, 8 



11 



3. Bond J.R., Cole S., Efstathiou G., Kaiser N., 
1991, ApJ, 379, 440 

4. Bullock J.S., et al., 2001, MNRAS, 321, 559 

5. Burigana C, Danese L., De Zotti G., 1991, 
A&A, 246, 59 

6. Burigana C, Finelli F., Salvaterra R., Popa 
L.A., Mandolesi N., 2004, in Recent Research 
Developments in Astronomy & Astrophysics 
- Vol. 2, 59 

7. Burigana C, Salvaterra R., 2003, MNRAS, 
342, 543 

8. Danese L., De Zotti G., 1980, A&A, 84, 364 

9. Dawson K.S., et al., 2002, ApJ, 581, 86 

10. De Zotti G., et al, 2004, in Proc. Int. 
Symp. Plasmas in the Laboratory and in the 
Universe: new insights and new challenges, 
Como, September 2003, eds. Bertin G., Fa- 
rina D., PozzoliR., pg. 375, |astro-ph/0401191 

11. da Silva A.C., et al., 2001, ApJ, 561, L15 

12. Ferrarese L., 2002, in Current high-energy 
emission around black holes, eds. C.-H. Lee 
& H.-Y. Chang, Singapore: World Scientific 
Pub. 

13. Fixsen D.J., et al., 1996, ApJ, 473, 576 

14. D.J. Fixsen and J.C. Mather, 2002, ApJ, 581, 
817 

15. Gnedin N.Y., Jaffe A.H., 2001, ApJ, 551, 3 

16. Granato G.L., et al., 2001, MNRAS, 324, 757 

17. Haiman Z., Loeb A., 1998, ApJ, 503, 505 

18. Hummer D.G., Storey P.J., 1987, MNRAS, 
224, 801 

19. Ikeuchi S., 1981, PASJ, 33, 211 

20. Jones M.E., 2003, in The scientific promise 
of the Square Kilometer Array, SKA 
Workshop, Oxford, 7 November 2002, 
eds. Kramer M., Rawlings S., pg. 47, 
http://www.skatelescope.org/documents/ 
Workshop- Oxford2002.pdf 

21. Lapi A., Cavaliere A., De Zotti G., 2003, ApJ, 
597, L93 

22. Knox L., 1995, Phys. Rev. D., 48, 3502 

23. Kogut A., 1996, Diffuse Microwave Emis- 
sion Survey, in XVI Moriond Astro- 
physics meeting Microwave Background 
Anisotropics, March 16-23, Les Arcs, France, 
astro-ph/9607100 

24. Kogut A., 2003, in the Proc. of The Cos- 
mic Microwave Background and its Polariza- 



tion, New Astronomy Reviews, eds. Hanany 
S., Olive K.A., Vol. 47, Issues 11-12, p. 945 

25. Kogut A., et al., 2003, ApJS, 148, 161 

26. Ma C.-P., Fry J.N., 2002, PRL, 88, 211301 

27. Mason B.S., et al., 2003, ApJ, 591, 540 

28. J.C. Mather et al, 1990, ApJ, 354, L37 

29. Mather J.C, Fixsen D.J., Shafer R.A., Mosier 
C, Wilkinson D.T., 1999, ApJ, 512, 511. 

30. Natarajan P., Sigurdsson S., Silk, J., 1998, 
MNRAS, 298, 577 

31. Natarajan P., Sigurdsson S., 1999, MNRAS, 
302, 288 

32. Navarro J.F., Frenk C.S., White S.D.M., 
1997, ApJ, 490, 393 

33. Oh S.P., 1999, ApJ, 527, 16 

34. Oh S.P., Cooray A., Kamionkowski M., 2003, 
MNRAS, 342, L20 

35. Ostriker J.P., Vishniac E.T., 1986, ApJ, 306, 
L51 

36. Partridge R.B., Richards E.A., Fomalont 
E.B., Kllerman K.I., Windhorst R., 1997 
ApJ, 483, 38 

37. Pei Y.,1995, ApJ, 438, 623 

38. Platania P., Burigana C, De Zotti G., Laz- 
zaro E., Bersanelli M., 2002, MNRAS, 337, 
242 

39. Prandoni I., et al., 2001, A&A, 365, 392 

40. Press W.H., Schechter P., 1974, ApJ, 187, 452 

41. Rees M.J., Ostriker J.P., 1977, MNRAS, 179, 
541 

42. Rephaeli Y., 1995, ARA&A, 33, 541 

43. Rybicki G.B., Lightman A.P., 1979, Radia- 
tive Processes in Astrophysics, Wiley, New 
York 

44. Salvaterra R., Burigana C, 2002, MNRAS, 
336, 592 

45. Springel V., White M., Hernquist L., 2001, 
ApJ, 549, 681 

46. Sasaki S., 1994, PASJ, 46, 427 

47. Silk J., 1968, ApJ, 151, 459 

48. Stebbins A., Silk J., 1986, ApJ, 169, 1 

49. Subrahmanyan R., Ekers R.D., 2002, in the 
XXVIIth General Assembly of the URSI, 
August 17-24, Maastricht, The Netherlands, 
|astro-phy 0209569 

50. Sunyaev R.A., Zeldovich Ya.B., 1970, 
Ap&SS, 7, 20 

51. Sunyaev R.A., Zeldovich Ya.B., 1972, Comm. 



12 



Astrophys. Space Phys., 4, 173 

52. Toffolatti L., et al., MNRAS, 1998, 297, 117 

53. Tremaine S., et al., 2002, ApJ, 574, 740 

54. Vishniac E.T., 1987, ApJ, 322, 597 

55. White S.D.M., Rees M.J., 1978, MNRAS, 
183, 341 

56. Zeldovich Ya.B., Sunyaev R.A., 1969, 
Ap&SS, 4, 301 



